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ABSTRACT 

The Herbig Ae star HD 139614 is a group-Ib object, which featureless SED indicates disk flaring and a possible 
pre-transitional evolutionary stage. We present mid- and near-IR interferometric results collected with MIDI, 
AMBER and PIONIER with the aim of constraining the spatial structure of the 0.1-10 AU disk region and assess 
its possible multi-component structure. A two-component disk model composed of an optically thin 2-AU wide 
inner disk and an outer temperature-gradient disk starting at 5.6 AU reproduces well the observations. This is an 
additional argument to the idea that group-I HAeBe inner disks could be already in the disk-clearing transient 
stage. HD 139614 will become a prime target for mid-IR interferometric imaging with the second-generation 
instrument MATISSE of the VLTI. 
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1. INTRODUCTION 

Understanding the time evolution of circumstellar disks around low-mass and intermediate-mass stars is of the 
highest importance to improve our knowledge of (exo)planet formation. It is currently admitted that the pro- 
cesses transforming gas rich disks into mature debris disks with planetary systems involve grain growth and 
settlinhg, photo-evaporation by stellar UV photons or disk clearing by forming and gas-accreting planets. How- 
ever, the mechanisms by which planet embryos form in the disk and evolve towards massive gaseous Jupiters or 
small rocky planets is still not well understood. In order to progress in this direction it is important to detect 
and identify peculiar spatial structures in circumstellar disks such as rings, gaps, asymmetries that may reveal 
planetary signatures at a particular evolutionary stage. Coupled to the modeling of spectral energy distribu- 
tions, infrared long baseline interferometric observations is a unique tool to probe and test disk configurations 
at various ages. 

In this paper we present the results of the observing campaign on the young star HD 139614 and the subse- 
quent modeling work based on VLTI interferometric data taken with MIDI, AMBER and PIONIER. With the 
low-resolution (R~30) mid-infrared data, the warm dust emission in the 1-10 AU region of the circumstellar disk 
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Figure 1. Spectral energy distribution of HD 139614 showing the near-infrared excess and the rising mid-infrared flux 
characteristic of pre-transitional disks. These two features can be modeled with simple hot (T >1000 K) and warm 
(T ~200-400K) blackbody components. 


is probed, while the H and K near-infrared data are sensitive to the hotter dust closer to the star and down 
to the dust sublimation radius. Our modeling has been based in parallel on the implementation of geometrical 
disk models and radiative transfer calculations. For a detailed review on the physics of disks and their spatial 
structure, we direct the reader to topical literature.^’ ^ 

In Sect. 2 we present the context of our study and the object HD 139614. In Sect. 3 and 4 the observations 
are presented together with the qualitative results. In Sect. 5 the modeling work is introduced and its outcome 
is discussed in Sect. 6. 

2. HD 139614 SYSTEM IN THE CONTEXT OF GROUP-I AND GROUP-H DISKS 


Teff 

Av 

M/Mq 

log(L/L0) 

R/Rq 

Age (Myr) 

Sp. Typ 

d(pc) 

7600T300 

0.5 

1.76T0.2 

l.lOiO.15 

2.06T0.4 


A7V 

142±27 


Table 1. Stellar parameter of HD 139614 taken from Alecian et al. (2013)® . 


HD 139614 is an intermediate-mass Herbig star with well characterized properties. In Table 1, we report the stellar 
parameters as derived in Alecian et al. (2013)® . It is a relatively evolved pre-main sequence star associated to the 
Lupus/Ophiuchus cloud. The stellar energy distribution (SED) exhibits a significant excess to the photospheric 
flux from the near-infrared to millimeter wavelengths. Its group-Ib featureless SED"^ suggests significant disk 
flaring as well as possible pre-transitional characteristics identified as a near-IR excess accompanied by a dip at 
6 microns in the SED, then followed by a rising mid-infrared part (cf. Fig. 1). 

In the Herbig HAeBe classification of Meeus et al. (2001)^ , HD 139614 is tagged as a group Ib object, meaning 
that its SED can be decomposed into a power law and a blackbody components and that it presents no silicates 
emission feature around 10 /im. This class is opposed to the group H(b), in which the SED is described solely by 
a single power law. The index a or b refer to, respectively, the presence or absence of a significant 10 /jm silicate 
feature. 

This classification is interesting from an evolutionary point of view in which a link between the two configurations 
is proposed: group I sources would host young hydrostatic disks with flaring geometry that produce the large 
mid- and far-infrared flux. Then these objects would evolve towards group II sources, which exhibit a flatter 
outer disk due to the settling of larger coagulated dust grains into the mid-plane. This results into a bluer and 
steeper infrared SED. The role of the puffed-up inner rim® in shadowing part of the disk of group II sources is 
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Figure 2. UV coverage obtained for HD 139614 with PIONIER (left panel), AMBER (center panel) and MIDI (right 
panel). Note that only the coverage of the quadruplet Al-Gl-KO-11 is displayed. 


also advanced. 

Recent works at high angular resolution have evidenced the presence of large gaps in disks classified as group I 
objects® . In such gapped disks, significant amount of mid- to far-IR emission in the SED can arise from large 
vertical walls at the inner edge of the outer disks that are directly irradiated by the central star. In addition, the 
presence or lack of silicate emission feature in the SED can also be interpreted in terms of gap size within the 
region where small and warm silicate grains (~300K) in a region where they are efficiently emitting. In contrast, 
for group II objects with flatter disks no gap has been reported so far by spatially resolved observations, while 
none of the studied group II objects show any absence of the silicate emission feature. The difference in the 
spatial structure of the inner disk regions between group I and group II sources is still an open question, and 
it is beyond the goal of this paper to discuss in details the conditions of evolution of HAeBe’s circumstellar 
disk. It is however interesting to note how recent spatially resolved observations have triggered new hypothesis 
on the possible evolutionary connection between group I and group II sources. It is suggested that these two 
classes could result from a bifurcation in the evolution of a common flared and gapless disk progenitor® . As 
the formation of gaps in protoplanetary disks has notable implications on the process of planet formation, it is 
essential to gather additional knowledge on a larger sample of the so-called pre-transitional disks 

3. OBSERVATIONS AND DATA REDUCTION 

HD 139614 was observed in several campaigns with the VLTI from 2010 to 2013. The MIDR observations were 
carried out with the 8-m unit telescopes in visitor mode on April 26*^ 2010 and in service mode on April 14*^ 
and 18*^ 2011. The observations were conducted using the popular HIGH-SENS mode of MIDI, which offers a 
resolution of R~30 at 10 )im. We used the baselines UT3-UT4 and UT2-UT3 for these runs. In the near-infrared, 
we first conducted AMBER® observations in service mode on May 9*/i 2012 using the H-l-K mode in low spectral 
resolution. We used the UT1-UT3-UT4 triplet under good seeing conditions (~0.8"). Additional H band data 
have been acquired with PIONIER® using the 1.8-m auxiliary telescopes between March 2012 and July 2013 
to improve the UV coverage for our object. PIONIER can recombine quadruplets instead of triplets as in the 
case of AMBER. In addition to the squared visibilities V®, both AMBER and PIONIER deliver closure phases 
using the different available triplets. However the closure phase measurements are not directly exploited in this 
paper. All the observations were obtained at low airmass (~1.1). In Table 2 the logs of the different observa- 
tions are presented together with the corresponding calibrator stars. The data reduction was achieved using the 
different pipelines publicly available such as EWS/MIA^® for MIDI or the amdlib package release 3.0.5 provided 
by the JMMC*. The MIDI calibrated visibilities and correlated fluxes span the 8-13 /^m spectral range, while 
the low spectral resolution AMBER squared visibilities span the 2.06-2.46 /im range in 13 spectral channels. 
With PIONIER we obtain three visibilities per baseline in three spectral channel within the H band. Both mid- 
and near-infrared observations present a good enough UV coverage to constrain different position angles across 
HD 139614’s disk (see Fig. 2). 
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Date 

UT 

Baseline 

Calibrator 

Seeing(") 



MIDI 





EOV of 300 mas (40 AU) 



26/04/2010 

03:16:05 

UT3-UT4 

HIP 60979 
HIP 63066 

1.2 

26/04/2010 

04:02:23 

UT3-UT4 

HIP 76552 
HIP 76397 

1.1 

26/04/2010 

05:35:51 

UT3-UT4 

HIP 76397 
HIP 76552 

1.3 

14/04/2011 

07:58:08 

UT1-UT2 

HIP 76397 

1.4 

18/04/2011 

03:35:05 

UT1-UT2 

HIP 72010 
HIP 74395 

0.7 

18/04/2011 

05:15:58 

UT1-UT2 

HIP 76552 

1.3 



AMBER 





EOV of 60 mas (9 AU) 



09/05/2012 

05:27:00 

UT1-UT3-UT4 

HD 141702 

1.01 

09/05/2012 

06:27:00 

UT1-UT3-UT4 

HD 140785 

0.65 



PIONIER 





EOV of 200 mas (25 AU) 



25/03/2012 

08:55:41 

Al-Gl-KO-Il 


1.3 

06/06/2013 

01:27:00 

A1-G1-J3-K0 


1.1 

16/06/2013 

23:31:00 

DO-Gl-HO-Il 


1.0 

03/07/2013 

23:58:00 

A1-B2-C1-D0 


1.1 


Table 2. Log of the observations of HD 139614 with MIDI, AMBER and PIONIER. 


For purpose of modeling, we also compiled the SED of our objects by retrieving optical photometry from the 
Tycho-2 catalog, near-infrared photometry from 2-MASS^^ , the SPITZER/IRS 5-40 /tm spectrum, the Her- 
schel/PACS photometry at 70, 100 and 160 /im, and millimeter photometry at 0.8 mm, 1.1mm, 1.3 mm and 
2 mm from Sylvester et al. (1997)^^ . 

4. OBSERVATIONAL RESULTS 
4.1. Spectral Energy Distribution 

The SED of HD 139614 is essentially featureless and presents no clear indication of silicate emission around 
10 /iin as inferred from the MIDI and the SPITZER/IRS spectra (see Eig. 3). The ISO spectrum published by 
Meeus et al. (2001) is displayed for comparison. The calibrated MIDI spectrum is a by-product of the EWS data 



Figure 3. Comparison between the MIDI, Spitzer and ISO spectra^ of HD 139614. The MIDI spectrum agrees well with 
Spitzer’s one and exhibit the PAH emission feature at 11.3 fira. 


reduction, and the error on the flux measurement was reduced by averaging the six flux measurements obtained 
for each UV point. The MIDI and Spitzer spectra agree very well over the whole N band, which suggest that 
no significant flux variability is detected. The PAH emission features seen the Spitzer spectrum at 8.6 /im and 
11.3 /rm are detected in the MIDI spectrum as well, although with less accuracy. In particular, the 11.3 /rm peak 
is clearly visible in the MIDI spectrum (see Sect. 6). 

4.2. Visibilities 

The spatial resolution achieved with our MIDI mid-infrared interferometric observations corresponds to ~20 mas 
(~3AU) at 10 /rm. In the near-IR, the observations are obtained both with the UTs (AMBER) and the ATs 
(PIONIER) baselines, spanning resolutions from ~2 to 11 mas (0.3-1.6AU). The calibrated visibility curves of 
MIDI shown in Fig. 4 indicate that the source is fully resolved for all the UV points with visibilities varying 
between 0.05 and 0.2 that show sinusoidal variations. As the hypothesis of a resolved close-binary system is 
discarded^^ such low visibilities indicate a significantly extended warm emitting region around HD 139614. 




Figure 4. Calibrated MIDI visibilities obtained with the different baselines shown in Fig. 2 (in black with error bars). 
The uniform achromatic ring model is plotted in orange, the one-component disk model in blue, and the two-component 
disk model in red for the high-density case (see Sect. 5.1). 


Furthermore the curves do not present the characteristic drop-off between 8 and 9 /im as seen for several Herbig 
stars disks in Leinert et al. (2004)^^ and modeled by Van Boekel et al. (2005)^® as the mid-infrared tail of the 
confined puffed-up inner rim, which implies high visibilities at 8 /rm followed by a sharp drop down to 9.5 /im. 
The emission at 8 /rm appears quite resolved by MIDI, suggesting that if hot material heated by the star is 
present it must be relatively extended. In the near-infrared the UT baselines mainly probe the high spatial 
frequencies and the hot dust in the system. As it can be seen in Fig. 5, the AMBER squared visibilities are quite 
low as well, i.e. varying between ~0.1 and 0.4, which suggest hot material extending well beyond the 0.3 AU 
distance theoretically resolved by the instrument. The PIONIER data also show resolved emission even with the 
shorter baselines (<40m) of the AT array, confirming that the hot dust emission is probably spatially extended. 
Although not shown in this article, both the MIDI differential phases and the AMBER and PIONIER closure 
phases did not reveal any clear spatial asymmetry in the spatial structure of the disk. 


5. MODELING 

5.1. MIDI and spectral energy distribution data 

The adopted methodology consisted in modeling the structure of the circumstellar disk from the spectral energy 
distribution and the VLTI interferometric visibilities. This is a commonly adopted approach that has demon- 
strated to be performant. We proceeded gradually from the modeling of the large scale parameters of the disk 
(inclination, P.A. etc...) using geometrical flat models to the more sophisticated 3D modeling of the inner disk 
using the RADMC-3D radiative transfer code^® . As a first step, we constrained the spatial extent of the 10 fira 
warm emitting region using solely the mid-IR visibilities. We applied a temperature-independent uniform ring 
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Figure 5. AMBER calibrated squared visibilities as a function of the wavelength for the six UV points (left) and PIONIER 
visibilities as a function of baseline at 1.6 /im (right). 


model parametrized by the inner radius of the ring its width Ai?, inclination i and position angle 0, and the 
photospheric contribution fin to the total flux^^ . The best fit obtained with two different grid resolutions^^ 
suggests that the bulk of the mid-infrared emission originates from a region extending from 4.5AU±0.2AU to 
8AU±0.2AU. A low inclination of the ring of i=20°±2.1° is derived, in agreement with previously published 
measurements^®’^® . A position angle of 6=112°±9.3° is found. 

We then adopted a more physical temperature-gradient disk model characterized by the temperature and surface 
density distributions Tr = Ti„{r / and = Tiin{r / rin)~^ ■ The first scenario implied a one-component disk 

model (i.e. central star+temperature-gradient disk), which parameters were optimized to fit simultaneously the 
MIDI visibilities and the near- and mid-infrared spectral energy distribution. The stellar SED was modeled with 
a Kurucz synthetic spectrum with Teff=7750K. A number of parameters were fixed in advance. These are the 
power law exponent of the surface density that was assumed to be p=3/2, a pure iron-free olivine composition^® 
and the outer radius fixed to 7?out=100 AU. For the disk density Ei„ at Hn we assumed three different cases of 
low (2xl0“®g.cm“^), medium (fix 10“®g.cm“^) and high (lO^'^g.cm"^) density to span a range of optical depth 
values (0.01 <T< 0.1) that allows the overall shape of the mid-infrared spectrum to be reproduced. We assume 
i=20° and 0=112° from the previous analysis. The free parameters for the one-component model are thus the 
inner radius rin, the temperature Tin at nn, and the temperature law exponent q. A simultaneous modeling of 
the SED and MIDI visibilities with x^ minimization resulted in the models shown in Table 3. We see that while 
the position of the inner rim of the outer disk is well constrained around 2.8 AU, the corresponding temperatures 
are quite scattered and high for a distance of ~3 AU. Note however that for the binary star V892Tau, Monnier 
et al. (2008)^^ derived a temperature of T ~450K at the inner rim of the circumbinary disk, which was found 
to be as far as ~17AU from the central star. However the modeling of the authors assumed for the equivalent 
luminosity of 400 Lq for two B8V stars. In all three cases, the best-fit models fail to reproduce both the SED 
and visibilities, in particular because almost no near-infrared excess is predicted. Smaller values for r\n obviously 
increases the near-infrared contribution in the model but also strongly disrupts the visibilities fit by predicting 


Model 

Surface density 

Arhot (AU) 

log(e) 

Tin (AU) 

TinK 

q 

Xr 


low 

- 


2.79±0.0fi 

fifi9±24 

0.74±0.0fi 

3.2 

One-component disk 

medium 

- 


2.78A0.09 

478±12 

O.fifiiO.Ofi 

3.0 


high 

- 


2.77±0.09 

421±9 

0.62±0.0fi 

3.4 


low 

1.6fi±0.12 

-2.78±0.1 

S.fiOiO.lO 

445±33 

0.78±0.18 

l.fi 

Two-component disk 

medium 

2.04A0.21 

-2.87±0.12 

5.fil±0.12 

3fi5±19 

0.97±0.21 

1.04 


high 

2.13A0.12 

-2.87±0.0fi 

5.fi3±0.12 

332±9 

1.04±0.0fi 

0.95 


Table 3. Results of the disk modeling after x^ minimization based on the MIDI visibilities and near-|-mid-infrared spectral 
energy distribution with associated 3 -(t error bars. 


10 ‘"r 

SEO 

Central star 

One-component model 

IWo-component model 



Figure 6. Near- and mid-infrared SED of HD 139614 and corresponding fit for the one-component and two-component 
models presented in this study for the high-density case 

high values at 8 /rm resulting from unresolved hot material close to the star. As the uniform ring model already 
predicted, the outer disk inner rim could be located further away, between 4.5 AU and 8 AU. This suggests a 
more sophisticated spatial structure, with for instance the presence of a hot inner component responsible for the 
near-infrared excess. We hence tested a second temperature-gradient model by introducing an inner component 
spatially separated from the outer disk. Such a structure is found in gapped pre-transitional disks^^’ . 

For our two-component temperature gradient model, we assumed a close-in hot component with a constant 
temperature of T ~1500K that produces the 3-4 pm excess typically observed for Herbig stars^^ . The hot dust 
component is assumed to be as close as r=0.22 AU, the dust sublimation radius derived from the luminosity 
of HD 139614. In addition to the parameter already defined for the one-component model, we introduce two 
additional free parameters, namely the radial width of the hot dust Arhot and the emissivity factor c to include 
the effect of the dust optical depth. The results of the minimization are displayed in Table 3. In addition 
to an improved result for the reduced the agreement with the observed SED near-infrared is much improved 
(see Fig. 6). The fit with the MIDI visibilities is still excellent as it can be seen in Fig. 4. Our best solution based 
on the combination of the MIDI visibilities and SED data has two components: an optically thin (e 
inner component extending from the dust sublimation radius to about 2.3 AU; a near-infrared emission originat- 
ing from a spatially confined region is at the contrary ruled out by our modeling; an optically thick outer disk 
starting at ~5.6 AU with a temperature of 332 K. However, as in the one-component disk case, our temperature 
estimate still exceeds the upper limit temperature (270 K) expected at 5.6 AU, under the assumption of a pure 
silicate dust composition. 

5.2. Near-infrared observations with AMBER and PIONIER 

As the possible hot inner component predominantly contributes to the near-infrared flux, the MIDI observations 
are not ideally suited to constrain unambiguously the dust structure within 1-2 AU from the star. The AMBER 
and PIONIER observations in the K and H bands offer the opportunity to confirm and refine our previous results. 
Therefore, we have implemented a self-consistent radiative transfer modeling of the spatial structure of the disk 
components. The code RADMC-3D code was used in this work. Based on a Monte Carlo approach, the code 
computes the radiation field and the disk temperature distribution from a given dust density distribution, and 
provides SEDs and synthetic images at different wavelengths. The dust surface density E(r) and scale height 
H{r) are parameterized as E(r) = Eout(?’/?’out)*’ and H[r) = Hout(r/rout)^ , where Eout and Hout are the surface 
density and scale height at the disk outer radius rout, P and /3 the power-law exponent and flaring index. 

The properties of the outer disk (power law exponent p, disk mass of the outer component Afout^, inner radius 
of the outer disk Rj^, dust composition) are constrained by the far-IR and sub-mm SED data points. 


^Total gas+dust, assuming gas-to-dust mass ratio of 100. 




Figure 7. Top left: synthetic SED of HD 139614 obtained by radiative transfer with the parameters shown in Table 4. 
Top right: fit to the experimental mid-IR visibilities. Bottom: Images of the disk of HD 139614 as seen at 2 /im and 10 pra 
with the star removed. We adopted a linear scale. The inclination i=20° is taken into account here. 


Then assuming co-planarity with the outer disk, the radiative transfer modeling of the inner disk is performed 
using a minimization exploring a 10x10x10 grid on the three free parameters rout, inner disk total mass 
Mhot^ and surface density index p. This is bounded by different values of the fixed parameters as the grain 
size distribution, the dust composition, and the ratio (H/r)ont- The accurate position of the outer disks inner 
radius Rm and the corresponding scale height is determined in a third step based on the MIDI visibilities. The 
numerical results of our modeling are summarized in the Table 4. 


Inner hot component 

rout (AU) 

P 

Afhot/ Mq 

grain size dist. (/rm) 

composition 

(H/r)out 

2.57 (r,ub=0.22AU) 

0.64 i/3=2/7) 

8.8xl0-y 

5-20 

80% olivine, 
20% graphite 

0.1 at 2.6 AU 

Outer disk 

Rin (AU) 

P 

Mout/ Mq 

grain size dist. (/rm) 

composition 

(H/r)out 

6.3 (Rout=150AU) 

-1 (/3=l/7) 

1.3x10“^ 

0.1-3000 

100% olivine 

0.14 at 150 AU 


Table 4. Results of the radiative transfer modeling of the multi-component disk of HD 139614. 


In Fig. 7 we show the results of the radiative transfer modeling and compare it to the observational data. The 
synthetic spectral energy distribution matches very well with the photometry from the optical to the millimeter 
range. The interferometric data are also quite well reproduced in the mid-infrared (cf. Fig. 7, upper right), and 
in the near-infrared (cf. Fig. 5 with the AMBER and PIONIER data). Synthetic images at 2/im and 10 /rm of 
the disk around HD 139614 have been produced using RADMC-3D as shown in Fig. 7, where they are plotted 
in a linear scale and with the central star removed. Not seen in these images is the tenuous dust that may still 
survive in the gap and which is needed to reproduce the excess emission in the SED around 6 fim. An upper 
limit of Mdust,gap<10~^'^ Mq is derived, above which the modeled interferometric visibilities clearly diverge from 
the measured quantities. 


*With the same assumptions as for the outer disk 




6. DISCUSSION AND CONCLUSIONS 


We show in this paper that the combination of near-infrared and mid-infrared interferometric data can strongly 
constrain the spatial structnre of the disk aronnd HD 139614. In particular, the interferometric data have ruled 
out with good confidence the hypothesis of a continuous disk in hydrostatic equilibrium and inwardly truncated 
at the dust sublimation radius. This reaffirms the power of long-baseline interferometry as a prime tool for 
studying young disks and their evolutionary trends. 

• HD 139614 appears to host an additional example of pre-transitional disks, i.e. a gapped disk found 
around both T Tauri and Herbig stars that could originate from the clearing process triggered by a young 
protoplanet. The possibility that the gap is in reality a strongly shadowed area of the disk is remote, as 
this would imply a very large scale height of the near-infrared disk with a brightness impacting strongly 
the near-infrared visibilities and at the same time a substantial amount of cold dust radiating at millimeter 
wavelength. 

• Our simulations indicate than instead of a confined optically thick and narrow ring of dust, the near- 
infrared excess is better modeled by a extended and optically thin dust component of ~2.5 AU width with 
an increasing dust density outward to the disk (positive parameter p) in order to match the level of mid-IR 
visibilities around 8 pm. The near-infrared visibility values do not contradict this hypothesis. 

• Preliminary results indicate that the scale height of the outer disk at its inner edge plays a significant role 
in reprocessing the mid- infrared flux and that its temperature should be relatively high (~350K). Further 
work will better constrain this quantity. 

• While not discuss in this paper, the PAH emission feature at 11.3/im is detected in the MIDI uncorrelated 
spectrum. Although the level of uncertainty on the MIDI visibilities does not permit to draw firm conclu- 
sions, we have indications on the fact that the PAH emitting region should be comparable in size, or even 
a bit larger, than the mid-IR continuum emission region. 

• None of near- and mid-infrared interferometric closure phases exhibit noticeable signatures. It is very 
unlikely that HD 139614 could be in reality a tight binary system with a low-mass infrared companion. 

Our high-resolution study of HD139614, a star classified as a group-Ib HAeBe source, strongly suggests the 
pre-transitional nature of its circumstellar disk, with a low dust-density gap extending over few astronomical 
units, and a relatively spatially extended and optically thin hot inner component. It is interesting to place 
this result in the larger context of HAeBe and T- Tauri pre-transitional objects: one finds for instance that the 
gaps in HD139614 and HD100546 are surprisingly narrow in comparison to the T-Tauri counterparts. In the 
latter objects however, the inner hot component was determined to be spatially compact. A non-exhaustive 



Age (Myr) 

Spectral type 

n (AU)^ 

rout (AU) 

D (pc) 

Gap size (AU) 

HD 139614^^ 

8 

A7 

0.2 

2.6 

140 

3.7 

HD 10054625.26 

10 

B9 

0.26 

^4 

103 

9 

HD 1691425 

13.5 

A5V 

~0.1 

0.2 (?) 

145 

<23(?) 

IRS 48® 

15 

AO 

~0.1 

0.3 

120 

<63 

HD 97048® 

- 

B9.5 

0.3 

2.5 

158 

<34 

HD135344B® 

6 

F4V 

~0.1 

0.3 (?) 

140 

<30(?) 

T Cha27 

7 

K0-G8 

0.07 

0.11 

100 

12 

LkCa 1523 

3-5 

K5-K3 

0.015 

0.19 

140 

58 

Rox 4423 


K3 

0.25 

0.4 

120 

36 

UXTau23 

1 

G8 

0.15 

0.4 

120 

36 


Table 5. Preliminary compilation of some known pre-transitional objects already published. The ages are given here 
without error bars, but an uncertainty of about several Myr should be accounted for. Hence we only have a trend on the 
age of these objects the ones respect to the others. (1): in most cased n is identified with the dust sublimation radius. 


compilation of currently known pre-transitional disks is presented in Table 5. Note that the ages are indicative 
and show, accordingly to the literature, uncertainties of several Myrs. 

Because of the importance of gapped disked in the understanding of planet formation, our result brings additional 
insight into this debate and propose an additional candidate to better understand the origin of disk dispersal 
mechanisms in different classes of objects. HD 139614 will become a prime target for mid-IR interferometric 
imaging with the second-generation instrument MATISSE of the VLTI. 
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